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Photoproduction of heavy quarks and exclusive production of vector mesons in ultraperipheral 
proton-nucleus and nucleus-nucleus collisions depend significantly on nuclear gluon distributions. 
In the present study we investigate quantitatively the extent of the applicability of these processes 
at the Large Hadron Collider (LHC) in constraining the shadowing component of nuclear gluon 
modifications. 
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I. INTRODUCTION 

Ultraperipheral relativistic heavy ion collisions can ex- 
plore several aspects of particle and nuclear physics and 
have been extensively discussed in the literature (for a 
small sample of references, see e.g. [H4l3|). In a re- 
cent publication we explored lead-lead peripheral colli- 
sions at energies available at the CERN Large Hadron 
Collider to probe gluon distributions in nuclei Us- 
ing recent gluon distributions from global fits to data, we 
investigated the sensitivity of direct photoproduction of 
heavy quarks and exclusive production of vector mesons 
to varying strength of gluon modifications. This idea, 
originally proposed in Ref. Q, can be used to constrain 
nuclear gluon distribution from data on production of 
heavy quarks and vector mesons. 

In this article we extend our previous work by consid- 
ering both direct and resolved processes in the photopro- 
duction of heavy quarks in ultraperipheral proton-lead 
(pPb) and lead-lead (PbPb) collisions at the LHC. We 
also consider, as an addition to the previous treatment, 
exclusive production of vector mesons in pPb collisions. 
Our main goal is to investigate quantitatively the extent 
of applicability of these processes in constraining nuclear 
gluon modifications. In order to present a self-contained 
report we include all relevant results and discussions from 
our previous study. 

The paper is organized as follows: below we briefly 
discuss the photon flux generated by one of the nuclei 
(or proton), making reference to several previous pub- 
lications where it has been discussed in details. Sec. |TT] 
treats the photoprocesses considered in this study, as well 
as the input nuclear and photon parton distributions. In 
Sec. IIIII we present the results of our calculations and a 
brief comment on theoretical errors. Our conclusion is 
contained in Sec. IIVI 

For a given impact parameter b, the flux of virtual pho- 
tons with photon energy fc, d 3 N 1 (k, h)/dkd 2 b, is strongly 
dependent on the Lorentz factor 7. At the Large Hadron 
Collider (LHC) at CERN the Lorentz factor in the labo- 
ratory frame 7l is 7455 for proton- proton (pp), 4690 for 
proton-lead (pPb) and 2930 for lead-lead (PbPb) colli- 
sions. The relationship between the Lorentz contraction 



factor associated with the relative velocity between the 
colliding nuclei, and the collider energy per nucleon, E / A, 
in GeV, is given by 7 = 2 7 £ - 1 w 2(1.0735£/^) 2 . The 
photon flux also de_p_e_nds strongly on the adiabaticity pa- 
rameter £ = kb/j 
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which drops off exponentially for £ > 1, above a cutoff 
energy determined essentially by the size of the nucleus, 
Ecutoff ~ 7 MeV/6 (fin). 

For symmetric nucleus-nucleus (AA) collisions, such 
as PbPb collisions at the LHC, each nucleus can act 
equally as source or target of the photon flux. Inte- 
grating d 3 N 7 (k,h)/dkd 2 b over impact parameters with 
the constraint of no hadronic interactions and account- 
ing for the photon polarization yields the total photon 
flux dN^(k)/dk given by 0, H3, 
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with Ra the radius of the nucleus. 

In the case of proton-nucleus (pA) collisions, the nu- 
cleus acts preferentially as the source and the proton as 
the target, leading predominantly to 7 p processes. But 
there is still a non-negligible contribution from 7A pro- 
cesses in which the proton acts as the source of photons 
and the nucleus as the target. Thus expressions for both 
types of fluxes are required for pA collisions. The flux 
due to the nucleus (of charge Z) can be evaluated ana- 
lytically and is given by [2J, 
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with reduced adiabaticity parameter, ( R , given by 
k(R p + Ra)/j and R p the effective radius of the 



proton. 
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The flux due to the proton is usually estimated using 
the dipole formula for the electric form factor [l5j : 
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where D = 1 + [0.71 GeV /Q^ in ] and the minimum mo- 
mentum transferred Q^m — ^ 2 /[7 2 (l — 2fc/V<SVjv)]- 

With the knowledge of the photon flux, any generic 
total photoproduction cross section can be expressed as 
a convolution of a process cross section u^(k) and the 
photon flux, dN^/dk. Thus for AA collisions 



Here H stands for a proton or a nucleus {H = p, A) and 
the total photoproduction cross section is obtained by 
convoluting the equivalent photon flux, dN 1 (k) / dk, with 
al H^QQX ^ ag described by Eq. © and Eq. (JBJ for 
A A and pA collisions. 

Let us now consider both processes in some detail. At 
leading order (LO), photon-gluon fusion leading to the 
production of a heavy quark pair is the only subprocess 
relevant to direct photoproduction. In view of the high 
energies involved, perturbative QCD is applicable, and 
the direct photoproduction cross section can be expressed 
as a convolution of the partonic cross section for the sub- 
process jg — > QQ and the relevant nucleon or nuclear 
gluon distribution: 
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with dN^ I dk given by Eq. ^ . The factor of 2 accounts 
for the source/target symmetry present in A A collisions. 
In the case of pA collisions, we have 
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with dN^ /dk and dN^/dk given by Eq. ([3]) and Eq. Q 
respectively. Relevant expressions for the rapidity distri- 
butions for both A A and pA collisions are given later. 



II. PHOTON-HADRON INTERACTIONS IN 
ULTRAPERIPHERAL COLLISIONS 



7 H->QQX 
direct 



(s) = |^a^W(s)/f(x,Q 2 )9(C), (8) 



with x the momentum fraction carried by the gluon. 
We use rriQ for the mass of the heavy quark (charm or 
bottom), s = W 2 H denotes the square of the center-of- 
mass energy of the photon-nucleus (or photon-nucleon) 
system, s — W 2 g that of the photon-gluon system, 
s — Am 2 -, 



and £ 



The gluon distribution in H, 



f^(x,Q 2 ), is evaluated at the pQCD factorization scale 
Q 2 = W 2 g = s. The function 9(C) enforces a minimum 
("threshold") value of x, x m i n , on the integral given by 
x mm = Am 2 q /W 2 H 

The partonic photon-gluon fusion cross section is given 
by [Ml 



A. Photoproduction of heavy quarks 

From the viewpoint of a Fock space decomposition, 
photon interactions with hadrons and nuclei can be clas- 
sified as direct or resolved. In direct interactions the 
photon behave as a point-like particle ("bare photon") 
while in resolved interactions the photon fluctuates into 
a quark-antiquark state or an even more complex par- 
tonic configuration consisting of quarks and gluons. The 
cross section for the photoproduction of a pair of heavy 
quarks, a^ H ^'QQ x (fc), is thus a sum of both the direct 
and resolved contributions, 
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with eQ the electric charge of the heavy quark Q, a. 



the electromagnetic coupling constant, (3 — Am^/s, and 



The resolved contribution is identical to the hadropro- 
duction of heavy quarks, and at leading order, involves 
terms corresponding to gluon-gluon and quark-antiquark 
subprocesses. The resolved cross section can be written 
as 
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where dx = dx±dx 2 , f^(xi, Q 2 ) (f^(x 2 , Q 2 )) is the dis- tribution of parton a with momentum fraction x\ (x 2 ) in 
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a photon (H) respectively, s = X1X2S and Q — s — 4mg. 
The summation over q runs over the light flavors, i.e. 
q — u,d, s. 

The partonic cross sections a"~*QQ (s) and a q ^®® (s) 
are given by [l9l - [2~J | 
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Here /3 and ^ are as defined previously. 

It is often more enlightening to represent the cross sec- 
tion in terms of rapidity, that is, to consider rapidity 
distributions. The differential cross section with respect 
to the rapidity of the heavy quark pair, da/dy, is related 
to the differential cross section with respect to photon 
energy, da/dk, through the relation da/dy = kda/dk. 
The rapidity distribution can therefore be expressed as 
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and scales directly with the photon flux dN^/dk. Thus 
for pA collisions, with the convention that the proton is 
incident from the right and the nucleus from the left, the 
total rapidity distribution is 
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where ki [ki cx e v ) and k r (k r cx e y ) simply denote 
photons from the nucleus and proton respectively, and 
the fluxes dN^ /dk and dN^/dk are given by Eq. ((3]) 
and Eq. (H|). As remarked earlier, the fluxes have sup- 
port only at small values of k, dying out exponentially 
at large k. Thus the first term on the right-hand side 
(7P distribution) peaks at positive rapidities while the 
second term (7A distribution) peaks at negative rapidi- 
ties. Since both the fluxes and process cross sections are 
different, the total distribution is manifestly asymmet- 
ric, and the 7p term dominates due to the much larger 
nuclear flux dN^ /dk. 

The total rapidity distribution for AA collisions can 
likewise be written as 



da AA^QQX 

dy 
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with ki (k r ) simply denoting photons from the nucleus in- 
cident from the left (right) and the flux dN^ / dk given by 
Eq. ([2]). Here the process cross sections and the left /right 
fluxes are identical, thus the respective rapidity distribu- 
tions are mirror images of each other, and consequently 
the total distribution is symmetric about midrapidity 
(0 = 0). 



B. Exclusive production of vector mesons 

The cross section for the exclusive elastic photoproduc- 
tion of a vector meson V on H (H = p, A) can be written 

as 



a 



lH^VH {k) _ 
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dt\F H (tf 
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where da 1 / dt\t—Q is the forward scattering ampli- 
tude and Fu(t) is the form factor. The dynamical infor- 
mation is encoded in the forward scattering amplitude 
while the momentum transfer of the elastic scattering is 
determined by the form factor, which is, in general, de- 
pendent on the spatial attributes of the target H . 

There have been studies of the photoproduction of J/\l/ 
and T in ultra per i p heral collisions at LHC (see for in- 
stance i, 0, [Hill I2H2I ) • In this work we use the sim- 
ple amplitude calculated from leading order two-gluon 
exchange in perturbative QCD [29|, |30| and corrected 
for other relevant effects (such as relativistic corrections, 
inclusion of the real part of the scattering amplitude, 
next-to-leading order NLO effects, etc, see for instance 
[3l], HH) through a phenomenological multiplicative cor- 
rection factor Cy- 

For elastic photoproduction on protons, the corrected 
LO scattering amplitude can be written as 



da lP^Vp 



dt 



Cv 



t=o 



16ir 3 a 2 s T e 
3aM v 



[xg p (x,Q 2 )] 2 . (17) 



Here, My is the mass of the vector meson (J/vP and 
T(ls) in the present study), x = M v /W 2 p is the frac- 
tion of the nucleon momentum carried by the gluons, 
and g p (x, Q 2 ) is the gluon distribution in a proton, eval- 
uated at a momentum transfer Q 2 = (My/2) 2 . T ee is 
the leptonic decay width and a s (a) the strong (electro- 
magnetic) coupling constant. Eq. [T7]is easily generalized 
to the nuclear case, 



da^ A ~ 
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[xg A (x,Q 2 )] 2 , (18) 



k=k r 



where gA(x,Q 2 ) = g p (x,Q 2 ) x R A (x,Q 2 ) is the nuclear 
gluon distribution and Rf(x, Q 2 ) the gluon modification. 

The correction factor Cy is estimated by constraining 
the calculated cross sections for elastic vector meson pho- 
toproduction on protons, a~ fp ^ Vp (W lp ), to reasonably 
reproduce the photoproduction data from HERA: [33[ 
for J/ty and |33l - l35| for T(ls). Further details can be 
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found in 14j. In contrast to photoproduction of heavy 
quarks, the quadratic dependence of the cross section on 
the gluon distribution has the significant implication of 
making exclusive vector meson production a very sensi- 
tive probe of nuclear gluon modifications. 

We now discuss the momentum-squared transferred 
dependence (^-dependence) of the cross section. For a 
proton the standard practice is to parametrize the t- 
dependence in the form of a rapidly decreasing exponen- 
tial function, e 6 '*', where b is referred to as the slope 
parameter. Integrating this function with respect to t 
gives a multiplicative factor 1/6 with units GeV 2 . In 
the case of a nucleus, the nuclear form factor FA(t) is 
given by the Fourier transform of the nuclear density 
distribution: FA(t) = J d 3 r Pa{t) e lq r , where q is the 
momentum transferred. For a heavy nucleus it is cus- 
tomary to model the density distribution as a Woods- 
Saxon distribution with parameters from electron scat- 
tering, pA(r) = Po/[l + e^ r ~' R - 4 ^ d l], with central density 
po, radius Ra and skin depth d. For 208 Pb in use at the 
LR(l p = 0.16/fm 3 , R A = 1.2A 1 / 3 fm, and d = 0.549 
fm |3. 

For the proton the photonucleon cross section 



(fc) is obtained through 



(T 



(k) = 



1 da^^P 



dt 



(19) 



with slope parameter b. In the present study we employ 
b = 4.5 GcV~ 2 . The photonuclear cross section is given 
by 



da"> A - 



■VA 



dt 



dt\F(t)\ 2 (20) 



t=0 Jt min (k) 



Here t n 



fm.m(k) = (M 2 /4fc7i) 2 , as is appropriate for narrow 
resonances [27} . The total cross sections for pA and AA 
collisions are obtained by a convolution of (j lp ^ Vp (k) and 
cf 1 (fc) with the relevant photon flux as described by 
Eq. ([6]) and Eq. ([5]) respectively. 

The photon energy, k, is related to the rapidity of the 
vector meson y by k = (My/2) exp(y). Thus, as in the 
case of photoproduction of heavy quarks, the differen- 
tial cross section with respect to rapidity is given by 
dal H^VH j Ay = (kdN 7 (k)/dk)a"t H ^ VH (k), and the cor- 
responding expressions for rapidity distributions in pA 
and AA collisions are of the forms given by Eq. (fT4")) and 
Eq. (|15p respectively. The symmetry attributes of the 
distributions are also similar. 



C. Parton distributions in nuclei and photons 

We now turn to the consideration of the parton dis- 
tributions (PDs) relevant to the processes considered in 
the present study. As mentioned earlier, the direct con- 
tribution to the photoproduction of heavy quarks is de- 
pendent solely on the gluon distributions in protons and 
nuclei, same as for the elastic photoproduction of vector 



mesons. The resolved contribution, on the other hand, 
requires the distributions of light quarks and antiquarks 
and gluons in photons, protons and nuclei. Thus in ad- 
dition to the usual requirement of nucleon and nuclear 
parton distributions (nPDs), there is also the need for 
the relatively poorly known parton distributions in pho- 
tons (7PDS), thereby increasing the level of the theoreti- 
cal uncertainties in the calculation of photoproduction of 
heavy quarks. Although the direct contribution is domi- 
nant, this dominance does not vitiate the need to have a 
good control on the resolved contribution, especially for 
bb production where the resolved components are quite 
sizable. 

Let us first discuss nuclear parton distributions. It is 
rather well-known that the distributions of partons (i.e. 
quarks and gluons) in nuclei are quite different from the 
distributions in free nucleons, that is, they are "modi- 
fied" by the complex, many-body effects in the nucleus. 
These nuclear effects are usually parametrized in terms 
of "nuclear modifications" R A (x,Q 2 ) which in general 
depend on the parton specie (a), the nucleus (A), mo- 
mentum fraction x and scale Q 2 . The nuclear effects can 
be categorized based on different intervals in a;. At small 
values of a; (x < 0.04), we have the phenomenon gener- 
ally referred to as shadowing. This is a depletion, in the 
sense that in this interval, the distribution of a parton a 
in the nucleus is smaller compared to the corresponding 
distribution in a free proton, i.e. R A < 1. Antishadow- 
ing, which is an enhancement (R A > 1), occurs in the 
range 0.04 <x < 0.3. Another depletion, the classic 
EMC effect pl], is present in the interval 0.3 < x < 0.8, 
while for x > 0.8, the Fermi motion region, we have an- 
other enhancement. It is important to note that although 
both shadowing and the EMC effect (antishadowing and 
Fermi motion) correspond to depletion (enhancement), 
the physical principles and mechanisms governing these 
phenomena are quite different. Further details can be 
found in (3^S2 With the knowledge of R A (x,Q 2 ), nu- 
clear parton distributions can be expressed as a convo- 
lution of free nucleon parton distributions and nuclear 
modifications, i.e. f A (x, Q 2 ) = f a (x, Q 2 ) <8> R A (x, Q 2 ). 

While the determination of quark and antiquark distri- 
butions in nucleons and nuclei is in general a nontrivial 
task, that of gluons is even more problematic. Gluons 
are electrically neutral, and thus their distributions can- 
not be extracted directly from Deeply Inelastic Scattering 
(DIS) and Drell-Yan (DY) processes which account for 
the major part of the data used in global fits. Their dis- 
tributions are in general inferred from sum rules and the 
Q 2 evolution of sea quarks distributions. The situation is 
even worse in the nuclear case: the available data is much 
less than for nucleons, and there is the added complica- 
tion of a mass dependence. It is therefore not unusual 
for nuclear gluon distributions from different global fits 
to differ significantly, especially in the magnitude of the 
various nuclear effects (shadowing, antishadowing, etc). 
This is especially obvious at low Q 2 (i.e. around their 
initial starting scales) since evolution to high Q 2 tends 
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to lessen the differences. Earlier global analyses [43-46] 
relied heavily on fixed-target nuclear deep-inelastic scat- 
tering (DIS) and Drell-Yan (DY) lepton-pair production. 
Incorporation of data on inclusive hadron production in 
deuteron-gold collisions has been implemented in [U |48[ , 



and neutrino-iron data in 
tion the approach in [52l . 



45 lj . We should also men- 
utilizing the Gribov pic- 
ture of shadowing. Despite all these advances the nuclear 
gluon distribution is still currently the least constrained 
aspect of global fits to nuclear parton distributions, as 
significant uncertainties still persist at both small and 
large x. 

Four recent nucleon and nuclear parton distributions 
are utilized in the present study. For the proton we 
use the Martin-Stirling-Thorne- Watts (MSTW08) par- 
ton distributions [54| which are available up to next-to- 
next-to- leading order (NNLO). In the nuclear case we 
use three nuclear modification sets. Two sets are by 
Eskola, Paukunnen, and Salgado, namely EPS08 and 
EPS09 [H Si|. The third is the Hirai-Kumano-Nagai 
(HKN07) distributions [H. While EPS08 is only to lead- 
ing order (LO), both EPS09 and HKN07 are available up 
to next-to- leading order (NLO). The distributions from 
MSTW08 serve two purposes: as the free nucleon distri- 
butions used in conjunction with nuclear modifications, 
and also as a "special" nuclear distribution in the absence 
of nuclear effects. The latter case is particularly useful 
for highlighting the influence of the various nuclear effects 
on observables. 

It is instructive to compare the characteristics of the 
gluon distributions from the four aforementioned sets 
based on the strength of their nuclear modifications. In 
Fig. [T] we show the nuclear modifications for gluons in Pb, 
flf b(x, Q 2 ) from EPS08, EPS09, and HKN07 at the fac- 



torization scale Q 2 = M 2 ^, appropriate for the elastic 
photoproduction of the J/\l/ meson. As already stated, 
one can view the MSTW08 gluon distribution as a nu- 
clear gluon distribution in the limit of zero nuclear effects 
( =1). At this scale, HKN07 has a rather weak gluon 
shadowing which extends well into the antishadowing re- 
gion, no antishadowing and gluon EMC effect, and an 
early onset of Fermi motion. EPS09 exhibits a moder- 
ately strong shadowing, and appreciable antishadowing 
and EMC effect, with a quite strong Fermi motion. Nu- 
clear modifications are strongest in EPS08: an especially 
strong shadowing, and substantial antishadowing, EMC, 
and Fermi motion. Thus in terms of shadowing we have 
a progression from zero effects to weak effects, moder- 
ate (intermediate) effects, then to strong effects as one 
progresses from MSTW08 to EPS08. The issue of un- 
certainties in parton distributions is becoming increas- 
ingly important especially in relation to precision tests 
of QCD. Discussions on uncertainties in nuclear parton 
distributions can be found in 0, |48| . 

Parton distributions in photons (7PDS) are derived 
from experimentally determined photon structure func- 
tion F^{x, Q 2 ), in conjunction with appreciable theoret- 
ical inputs. These inputs, which are necessary in im- 
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FIG. 1: (Color online) Nuclear gluon modifications in Pb, 
R™{x,Q 2 = Mj /9 ), from EPS08 (dashed line), EPS09 (solid 
line), and HKN07 (dash-dotted line) respectively. 



plementing the parametrizations of photon parton distri- 
butions from the structure function, account in part for 
some of the observable differences in the available pho- 
ton parton distribution sets. Another source of differ- 
ences is in the choice and scope of experimental data from 
which F2 is extracted. At present there is an apprecia- 
ble number of photon parton distribution sets available, 
both at leading and next-to-leading orders [55j-t70j. It 
should be noted that unlike in the case of a nucleon, there 
are no valence quarks present in the photon; therefore 
antiquark distributions are the same as quark distribu- 
tions. Furthermore, there are no sum rules governing the 
gluon content; thus the gluon distribution is almost to- 
tally unconstrained. The gluon distribution contributes 
to F2(x,Q 2 ) majorly through the j*g — > qq channel, 
which has significant numerical support only at small x. 

In light of the relatively poor situation especially in 
the gluon sector, we use three different photon parton 
distribution sets in the present study. The rationale for 
this will become clear when we present results for the 
resolved component in heavy quark production. Two 
sets are relatively older and have been in common use: 
the Gluck-Reya-Vogt (GRV) set [H| and that by Schuler 
and Sjostrand (SaSld) (6a |. The third set, the Cornet- 
Jankowski-Krawczyck (CJK2) set [6!|, is more recent and 
has facility to estimate error on a calculated quantity due 
to the uncertainties in the input photon parton distribu- 
tions. In Fig. [2] we show the parton distributions from 
these three sets for the light (u,d,s) quarks and gluon. 

The up quark is the best constrained component of 
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evaluated to one loop at the scale Q 2 using the evolution 
code contained in the MSTW08 package. Two comments 
pertaining to the results on photoproduction of heavy 
quarks are in order: first, due to the relative smallness 
of nuclear effects in heavy quark photoproduction, espe- 
cially for pPb collisions, we have, for visual purposes, left 
out the results from the HKN07 distributions. In gen- 
eral the HKN07 values are intermediate between those of 
EPS09 and MSTW08. Second, while we have used three 
sets of photon parton distributions in the calculations, 
we show only the rapidity distribution plots using the 
GRV set, since the GRV results are consistently some- 
what intermediate between those from CJK2 and SaSld 
distributions. 



A. Heavy Quarks in pA and AA Collisions 

1, Open charm 



FIG. 2: (Color online) Parton distributions in the photon for 
light quarks (u,d,s) and gluons at Q 2 = 10 GeV 2 for three 
distribution sets: CJK2 (solid lin e) [69| . GRV (dashed line) 
[Hi]], and SaSld (dash-dotted line) [65[ respectively. For better 
visuality the down quark distributions have been enhanced by 
a factor of 5/3, the strange quark distributions by 5/2, and 
the gluon distributions by 1/14. 



photon parton distributions, and even in this case there 
are clear differences between the three sets. This is also 
the case for the down quark. Only for the strange quark 
distribution is there good agreement between the three 
sets. But these differences pale in comparison with that 
of gluons where the differences between the sets are quite 
significant almost across the whole x interval, with CJK2 
having the largest magnitude while SaSld has the lowest, 
with the GRV somewhat intermediate. For the LO re- 
solved contribution to heavy quark photoproduction, the 
two competing partonic subprocesses are the gg —> QQ 
and qq — !• QQ, and it thus to be expected that signifi- 
cant differences will arise in the resolved cross sections 
whenever the gg — > QQ component is appreciable. 



III. RESULTS 

We now discuss the results of our calculations for both 
the inclusive photoproduction of heavy quarks (cc and 
bb) and the exclusive elastic production of vector mesons 
(J/* and T(ls)) in pPb ( V /5^T = 8.8 TeV) and PbPb 
{^/snn — 5.5 TeV) collisions at the LHC. We take 
m c = 1.4 GeV and nib — 4.75 GeV for consistency with 
the MSTW08 parton distributions, and mj^ = 3.097 
GeV and my = 9.46 GeV respectively. The strong cou- 
pling constant a s {Q 2 ), needed for the calculations, is 



In Table U we present the cross sections for both di- 
rect and resolved cc production in ultraperiphcral pPb 
collisions at the LHC, using the previously-described nu- 
cleon/nuclear and photon parton distributions. Due to 
the large disparity in the magnitude of the photon fluxes 
from the proton and the nucleus (Pb), it is a common 
practice to neglect the photonuclear (7Pb) contribution 
to the total cross section. Another simplification is the 
neglect of the resolved contribution which is relatively 
small in comparison with the dominant direct contribu- 
tion. Since both the 7Pb and resolved contributions are 
dependent on nuclear and photon parton distributions, 
the results presented in Table U afford a good handle in 
determining the degree of validity of the above approxi- 
mations for the rather significantly different attributes of 
the distributions utilized in the present study. 

TABLE I: Cross sections for photoproduction of cc in ultra- 
peripheral pPb collisions at the LHC. All cross sections are in 
microbarns (fJo). 





PDF 


Direct 


Resolved 


SaSld 


GRV 


CJK 


7P 


MSTW08 


5570 


692 


1157 


1418 




MSTW08 


607 


114 


195 


228 


7A 


EPS08 


376 


94 


162 


185 




EPS09 


471 


102 


175 


202 



We first consider the photonuclear part. The 7Pb 
contribution to the total cross section is approximately 
10.6% for MSTW08 (no nuclear modifications), 8.8% 
for EPS09 ( moderate gluon shadowing), and 7.4% for 
EPS08 (strong gluon shadowing), almost totally indepen- 
dent of the choice of photon parton distribution set. The 
resolved contribution, on the other hand, while depend- 
ing strongly on the choice of photon parton distributions, 
has an even more negligible dependence on the choice of 
nuclear parton distributions. For the SaSld distribution, 
the resolved component is approximately 12% while it is 
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approximately 18% for the GRV distribution, and 21% 
for the CJK2 set. As is to be expected about 96% of the 
resolved component stems from 7p while about 90% of 
the 7Pb contribution derives from the direct component. 
Overall the sensitivity to nuclear modifications (domi- 
nantly shadowing) is rather small: the no-modification 
(MSTW08) total cross section is reduced by approxi- 
mately 4% and 2% by the modifications in EPS08 and 
EPS09 respectively. 
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FIG. 3: (Color online) Rapidity distributions of cc photopro- 
duction in pPb collisions at the LHC using the GRV photon 
parton distributions. Left panel: 7p and 7Pb contributions to 
total rapidity distributions. Dotted line depicts the 7p con- 
tribution while the dashed (MSTW08), solid (EPS09), and 
dash-dotted (EPS08) lines correspond to 7Pb contributions 
with no shadowing, moderate shadowing, and strong shad- 
owing respectively. Right panel: total rapidity distributions 
(sum of 7p and 7Pb contributions). 

In Fig. [3] we show the corresponding rapidity distri- 
butions for cc production using the GRV photon parton 
distributions. The left panel depicts the distributions 
separately for both 7p and 7Pb contributions while the 
right panel shows the total distributions, i.e. the sum of 
both contributions. As is apparent from the left panel, 
the 7p contribution is dominant at all rapidities except 
for y < — 5. In line with our convention, the 7p contri- 
bution peaks at positive rapidities while the 7Pb distri- 
butions peak at negative rapidities, and the asymmetric 
nature of the total distributions is clearly exhibited. De- 
spite the smallness of the 7Pb contribution, the effects 
of nuclear modifications are still slightly manifested in 
the total rapidity distributions shown in the right panel. 
It is interesting to note that the sensitivity to nuclear 
effects is significantly enhanced when one considers the 



ratio of the total rapidity distribution at a specific neg- 
ative rapidity (— yi), and the rapidity distribution at 
the equivalent positive rapidity (+yi), somewhat rem- 
iniscent of pseudorapidity asymmetry in deuteron-gold 
collisions. The MSTW08 ratio for y = -3,3 is reduced 
by about 9% and 5% by the modifications in EPS08 and 
EPS09 respectively, while for y — —4, 4 the reductions 
amount to 15% and 10%, and 23% and 14% respectively 
for y = —5,5. Thus the rapidity asymmetry ratio in the 
interval — 5<j/<4,4<y<5 exhibits appreciable 
sensitivity to nuclear effects. 

We now discuss the corresponding case of cc produc- 
tion in ultraperipheral PbPb collisions, and in Table UU 
we show the cross sections (in mb) for both direct and 
resolved contributions. 

TABLE II: Cross sections for photoproduction of cc in ultra- 
peripheral PbPb collisions at the LHC. All cross sections are 
in millibarns (mb). 
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180 
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EPS08 


890 


99 


165 


207 


EPS09 


1002 


101 


169 


213 



The sensitivity to nuclear modifications (dominantly 
shadowing) is more significant here than in pPb colli- 
sions. Thus considering the direct contribution which 
is the dominant part, the modification- free (MSTW08) 
cross section is reduced by approximately 24% by the 
rather strong nuclear shadowing present in EPS08, and 
by about 14% in the case of EPS09. With the inclusion of 
the resolved component, the reductions are around 22% 
and 13% respectively, and vary only slightly with the 
choice of photon parton distribution. 

The resolved contribution exhibits a mild dependence 
on nuclear modifications and a much stronger sensitivity 
to photon parton distribution. Thus in the case of SaSld 
the resolved component is 8.6% of the total cross sec- 
tion for MSTW08, 10% for EPS08, and 9.2% for EPS09. 
For GRV, the contributions are 13.4%, 15.6%, and 14.4% 
respectively. The CJK2 values are the largest: 16.2%, 
18.9%, and 17.5% respectively. 

The rapidity distributions for cc production in ultra- 
peripheral PbPb collisions are shown in Fig. [H and are 
manifestly symmetric about midrapidity (y — 0). The 
sensitivity to nuclear modifications is more transparent 
here than in total cross sections. Shadowing is the dom- 
inant nuclear effect for — 3 < y < 3, and the rapidity dis- 
tributions in this region reproduce the observed trend of 
gluon shadowing strength exhibited in Fig. [TJ MSTW08 
with its zero gluon shadowing gives the largest rapidity 
distribution while EPS08, with its strong gluon shadow- 
ing, gives the smallest. Due to strong flux suppression, 
shadowing is most markedly apparent for the rapidity 
range ~2 < y < 2. This range therefore provides a good 
window to discriminate among different gluon shadowing 
scenarios. 
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FIG. 4: (Color online) Total rapidity distributions of the pho- 
toproduction of cc in PbPb collisions at the LHC using the 
GRV photon parton distributions. Dashed line depicts result 
using the MSTW08 gluon distribution (no nuclear modifica- 
tions). Solid and dash-dotted lines are results from nuclear- 
modified gluon distributions from EPS09 and EPS08 respec- 
tively. 



The rapidity intervals 3 < y < 6 corresponds to x m i n 
in the antishadowing region (deep shadowing) for right 
(left) incident photons and vice versa for —6 < y < —3. 
Due to the photon flux suppression in the deep shadowing 
region, the rapidity distributions are sensitive mainly to 
antishadowing in addition to both EMC effect and Fermi 
motion. Since both EPS08 and EPS09 have substantial 
antishadowing, their rapidity distributions reflect this, 
being slightly larger than those from MSTW08. The dis- 
criminatory power here is not as appreciable as in the 
shadowing case though, due largely to the smallness of 
the distributions. 

For both rapidity ranges y < —6 and y > 6, x m i n > 0.2 
and the relevant contributing nuclear effects are the EMC 
and Fermi motion since the contribution from antishad- 
owing is small, and that from shadowing practically 
nonexistent by virtue of flux suppression. Both EPS08 
and EPS09 nuclear modifications exhibit EMC effect and 
Fermi motion, and the destructive interference from both 
effects render their rapidity distributions to practically 
coincide with that from MSTW08. 



2. Open bottom 

We now discuss our results for total cross sections and 
rapidity distributions for bb production. In Table Mil we 



present the cross sections for both direct and resolved bb 
production in ultraperipheral pPb collisions at the LHC. 
In parallel with the treatment of cc production, we de- 
termine the relative importance of 7Pb and resolved con- 
tributions to the total bb production cross sections. For 



TABLE III: Cross sections for photoproduction of bb in ultra- 
peripheral pPb collisions at the LHC. All cross sections are in 
nanobarns (nb). 





PDF 


Direct 


Resolved 


SaSld 


GRV 


CJK 


7P 


MSTW08 


36512 


8641 


12178 


14977 




MSTW08 


5084 


2061 


3032 


3663 


7A 


EPS08 


3972 


1886 


2808 


3389 




EPS09 


4409 


1936 


2874 


3477 



MSTW08 (no nuclear modifications) the 7Pb contribu- 
tion to the total cross section is approximately 14%, 13% 
for EPS09 ( moderate gluon shadowing), and 12% for 
EPS08 (strong gluon shadowing), with little dependence 
on the choice of photon parton distribution set. It is thus 
apparent that while the photonuclear contribution is rel- 
atively more significant here than in cc production the 
effect of nuclear shadowing is less. The resolved contribu- 
tion again depends significantly on the choice of photon 
parton distribution, and exhibits negligible dependence 
on the choice of nuclear parton distribution. For the 
SaSld distribution, the resolved component is approxi- 
mately 21% while it is approximately 27% for the GRV 
distribution, and 30% for the CJK2 set. These values 
are larger than the corresponding ones for cc produc- 
tion, and thus serve to underscore the relatively higher 
importance of the resolved component in bb photopro- 
duction. As noted above, the sensitivity to nuclear mod- 
ifications (dominantly shadowing) is very small: the no- 
modification (MSTW08) total cross section is reduced 
by approximately 2.4% and 1.5% by the modifications 
in EPS08 and EPS09 respectively. Thus with its appre- 
ciably large resolved component and weak sensitivity to 
nuclear effects, it is tempting to speculate that bb photo- 
production could be of some use in constraining photon 
parton distributions. 

In Fig. [5] we show the corresponding rapidity distri- 
butions for bb production using the GRV photon parton 
distributions. The left panel depicts the distributions 
separately for both 7p and 7Pb contributions while the 
right panel shows the total distributions, being the sum 
of both contributions. As is apparent from the left panel, 
the 7P contribution is dominant at all rapidities except 
f° r V — 4. As in cc production the 7p contribution 
peaks at positive rapidities while the 7Pb distributions 
peak at negative rapidities, and the asymmetric nature 
of the total distributions is clearly manifested. Due to 
the smallness of the nuclear effects, the total rapidity 
distributions shown in the right panel almost overlap for 
the different nuclear parton distributions considered in 
the present study. Similar to cc production, the rapid- 
ity asymmetry ratio exhibits appreciable sensitivity to 
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FIG. 5: (Color online) Rapidity distributions of bb photopro- 
duction in pPb collisions at the LHC using the GRV photon 
parton distributions. Left panel: 7p and 7Pb contributions to 
total rapidity distributions. Dotted line depicts the 7p con- 
tribution while the dashed (MSTW08), solid (EPS09), and 
dash-dotted (EPS09) lines correspond to 7Pb contributions 
with no shadowing, moderate shadowing, and strong shad- 
owing respectively. Right panel: total rapidity distributions 
(sum of 7p and 7Pb contributions). 



component is approximately 15% of the total cross sec- 
tion for SaSld, 19% for GRV, and about 23% for CJK2. 
Again, these values indicate that the resolved component 
is relatively more significant in bb as compared to cc pho- 
toproduction. 

In Fig. [5] we show the rapidity distribution for bb in 
ultraperipheral PbPb collisions at the LHC. Shadowing 
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nuclear effects around y — —4, 4: the MSTW08 ratio is 
reduced by 19% and 6% respectively by the modifications 
in EPS08 and EPS09. 

Let us now consider the corresponding case of bb pro- 
duction in ultraperipheral PbPb collisions. In Table IIVI 
we show the cross sections (in /jb) for both direct and 
resolved contributions. 

TABLE IV: Cross sections for photoproduction of bb in ultra- 
peripheral PbPb collisions at the LHC. All cross sections are 
in microbarns (/ib). 
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The sensitivity to nuclear modifications, while more 
appreciable here than in pPb collisions, is signifi- 
cantly less than for the equivalent cc production. The 
modification-free (MSTW08) cross section is reduced 7% 
by the modifications in EPS08, and about 5% in the case 
of EPS09. The resolved contribution is almost indepen- 
dent of nuclear modifications but exhibits a strong sen- 
sitivity to photon parton distribution. Thus the resolved 



FIG. 6: (Color online) Total rapidity distributions of bb pho- 
toproduction in PbPb collisions at the LHC using the GRV 
photon parton distributions. Dashed line depicts result using 
the MSTW08 gluon distribution (no nuclear modifications). 
Solid and dash-dotted lines are results from nuclear-modified 
distributions from EPS09 and EPS08 respectively. 



dominates in the rapidity interval — 2 < y < 2 and is most 
clearly manifested in the rapidity window — 1 < y < 1. 
Thus this interval presents the best sensitivity to shad- 
owing effects in bb production. Although less marked, the 
progression of rapidity distribution with relative shadow- 
ing strength follows the trend observed in cc production: 
MSTW08 still gives the largest distribution while EPS08 
gives the smallest. As in the case of cc production, there 
is a slight manifestation of the influence of antishadowing 
around y = — 3 and y = 3. The distributions practically 
overlap for y < —4 and y > 4; thus overall, the interval 
— 1 i5 V ^ 1 seems to afford the best sensitivity to nu- 
clear effects, in this case primarily shadowing. A detailed 
treatment of bb production at the LHC with emphasis on 
the x dependence can be found in [71j . 
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B. Vector Mesons in pA and A A Collisions 

1. J/* 

We now present our results on elastic photoproduc- 
tion of the J/vP and T(ls) in the framework of a cor- 
rected leading-order two-gluon exchange formalism in 
QCD. Here the nuclear effect (shadowing) is quite siz- 
able; we therefore also display explicitly results from the 
HKN07 distribution set. 

In Table fVl we present the components and total cross 
sections for the elastic photoproduction of J/^S in ultra- 
peripheral pPb collisions at the LHC. As expected the 
7P contribution is dominant, with the 7Pb contribution 
decreasing in size with increasing severity of gluon shad- 
owing. 

TABLE V: Cross sections (in fib) for elastic photoproduction 
of J/\P in ultraperipheral pPb collisions at the LHC. Second 
and third columns are the contributions from 7p and 7Pb 
respectively for different distributions. The sums of the two 
contributions are presented in the fourth column. 





7P 7A Total 


MSTW08 
EPS08 
EPS09 
HKN07 


167.3 23.6 190.9 
2.2 169.5 
8.5 175.8 
15.4 182.7 



The contribution of the 7Pb component to the to- 
tal cross section ranges from a high of about 12% for 
MSTW08 (no gluon shadowing) to a low of 1.3% for 
EPS08 (strong gluon shadowing), in accordance with the 
trend displayed in Fig. [TJ This trend is replicated in the 
degree of shadowing reflected in the total cross section: 
the no-shadowing (MSTW08) cross section is reduced by 
about 11% by the shadowing in EPS08, 8% in EPS09, 
and 4% in HKN07 respectively. 

Fig. [7] shows the rapidity distributions for both 7p and 
7Pb components and their sums. The distributions are 
manifestly asymmetric, with the dominant 7p component 
peaking at positive rapidities while the subdominant 7Pb 
contributions peak at negative rapidities. The 7Pb dis- 
tributions exhibit clearly the influence of gluon shadow- 
ing, especially in the narrow rapidity window —3 < y < 1 
where the differences mimic the relative strength of gluon 
shadowing in the respective nuclear parton distribution. 
Although the 7p contribution is dominant for practically 
all rapidities, these 7Pb distributions with the excep- 
tion of EPS08 are appreciable enough in the said rapid- 
ity range such that the total rapidity distributions show 
the same tendency as that of the 7Pb component in the 
designated rapidity window. Thus it seems feasible that 
a consideration of J/^P production in pPb collisions in 
this rapidity interval offers some potential in constrain- 
ing gluon shadowing. 

We now turn to J/'J' production in utraperipheral 
PbPb collisions at the LHC and in Table I VII we present 
the total cross sections from the different parton distribu- 
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FIG. 7: (Color online) Rapidity distributions of exclusive pho- 
toproduction of J /\P in pPb collisions at the LHC. Left panel: 
7P and 7Pb contributions to total rapidity distributions. 
Dotted line depicts the 7p contribution while the dashed 
(MSTW08), dash-double- dotted (HKN07), solid (EPS09), 
and dash-dotted (EPS08) lines correspond to 7Pb contribu- 
tions with no shadowing, weak shadowing, moderate shadow- 
ing, and strong shadowing respectively. Right panel: total 
rapidity distributions (sum of 7p and 7Pb contributions). 



tions. The differences in the predicted cross sections are 
extremely clear cut, and thus this process can serve as 
an excellent probe of gluon shadowing as well as a good 
discriminator of the different gluon shadowing templates 
considered in the current study. 

TABLE VI: Total cross sections (in mb) for elastic photo- 
production of J/ty in ultraperipheral PbPb collisions at the 
LHC. 
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It is instructive to consider the remarkable effect that 
the quadratic gluon dependence has on the cross sections 
(and also on the rapidity distributions). Thus the no 
shadowing cross section is reduced by approximately 87% 
by the strong shadowing in EPS08, and by about 61% 
by the relatively moderate shadowing in EPS09. Even 
the rather weak gluon shadowing present in HKN07 is 
responsible for reducing the cross section by about 34%. 

The constraining potential of the J/*f? elastic produc- 
tion process is more clearly exhibited when rapidity dis- 
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tributions are considered. To this end we display the 
corresponding rapidity distributions in Fig.[5]for the four 
different gluon shadowing scenarios. 



E 
> 

TJ 
X3 







1 


1 1 1 1 1 


- MSTW08 






~-\ - EPS08 




/ 

/ 


\ - EPS09 




/ 

/ 

/ 

/ 


\ - HKN07 - 

\ 
\ 




/ 

1 / 

1 / 
/ i 

l: 


i 

\ i 
\ \ 

w 

M 




1! 

ji 






"/ 


\\\ 












V 









FIG. 8: (Color online) Total rapidity distributions of exclu- 
sive photoproduction of J/ty in PbPb collisions at the LHC. 
Dashed line depicts result using the MSTW08 gluon distribu- 
tion (no nuclear modifications). Solid, dash-dotted, and dash- 
double-dotted lines are results from nuclear-modified gluon 
distributions from EPS09, EPS08, and HKN07 parton distri- 
butions respectively. 

Shadowing is the relevant nuclear effect in the rapidity 
interval — 3 < y < 3 and unsurprisingly, the rapidity dis- 
tributions mimic the behavior in the shadowing region 
of Fig. [TJ The largest rapidity distribution is given by 
MSTW08, followed by HKN07, and EPS09. The smallest 
is by EPS08 due to its strong gluon shadowing. The ra- 
pidity window — 2 < y < 2 manifestly depicts the signif- 
icant distinction between the various gluon distributions 
arising from the quadratic dependence. Antishadowing 
manifests in the intervals —5 < y < — 4 and 4 < y < 5; 
the effect though is quite slight. 

2. T(ls) 

We now discuss our results for elastic production of 
T(ls) in ultraperipheral pPb and PbPb collisions. The 
ensuing treatment parallels closely that of J /'F in many 
respects due to the same underlying production mecha- 
nism. 

In Table Ivni we present the components and total cross 
sections for the elastic photoproduction of T in ultrape- 
ripheral pPb collisions at the FHC. While the 7p contri- 
bution is still dominant, the 7Pb contributions from the 



different gluon distributions are markedly significant, and 
exhibit the trend of decreasing magnitude with increasing 
strength of gluon shadowing. 

TABLE VII: Total cross sections (in nb) for elastic photopro- 
duction of T in ultraperipheral pPb collisions at the LHC. 
Second and third columns are the contributions from 7p and 
7Pb respectively for different distributions. The sums of the 
two contributions are presented in the fourth column. 





7P 


7A 


Total 


MSTW08 


390 


219 


609 


EPS08 




84 


474 


EPS09 




130 


520 


HKN07 




161 


551 



The contribution of the 7Pb component to the total 
cross section is 36% for MSTW08 (no gluon shadowing), 
29% for HKN07 (weak gluon shadowing), 25% for EPS09 
(moderate gluon shadowing), and 18% for EPS08 (strong 
gluon shadowing). This gradation accords with the trend 
observed in Fig. [TJ which is also replicated in the degree 
of shadowing reflected in the total cross section: the no- 
shadowing (MSTW08) cross section is reduced by about 
22% by the shadowing in EPS08, 15% in EPS09, and 
10% in HKN07 respectively. This appreciable magni- 
tude of the effect of shadowing seems to indicate that 
T photoproduction cross section in ultraperipheral pPb 
collisions offers some promising potential in constraining 
nuclear gluon shadowing. 

In Fig. [5] we show the rapidity distributions for both 
7P and 7Pb components (left panel )and the total (right 
panel). As usual the distributions are asymmetric, with 
the dominant 7P component peaking at positive rapidi- 
ties while the subdominant 7Pb contributions peak at 
negative rapidities according to the convention adopted 
in the present study. 

Let us consider the left panel. The 7Pb component 
distributions here, relative to that of 7p, are more appre- 
ciable than for J^> production. In fact, for — 5 < y < — 2 
these distributions are larger than the 7p distribution, 
and around their peaks at y ~ —2.7 they reflect quite 
distinctly the effect of the varying shadowing strength in 
the gluon distributions used. Thus in the rapidity inter- 
val —4 < y < — 1 the total distributions in the right panel 
show good sensitivity to gluon shadowing, and therefore 
afford good potential for constraining purposes. 

We now turn to T production in ultraperipheral PbPb 
collisions at the LHC. In Table rVTLTl we present the total 
cross sections from the different gluon distributions. 

Although not quite as dramatic as in the case of 
here also one can see the effect on the cross sections of 
the quadratic dependence on gluon distribution. The 
MSTW08 (no-shadowing) cross section is reduced by ap- 
proximately 48% by the strong shadowing in EPS08, 
by about 31% by the relatively moderate shadowing in 
EPS09, and by 23% by the rather weak gluon shadow- 
ing present in HKN07. Thus a consideration of the cross 
section can serve as a good probe of gluon shadowing as 
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FIG. 9: (Color online) Rapidity distributions of exclusive 
photoproduction of T(ls) in pPb collisions at the LHC. 
Left panel: 7p and 7Pb contributions to total rapidity dis- 
tributions. Dotted line depicts the 7p contribution while 
the dashed (MSTW08), dash-double- dotted (HKN07), solid 
(EPS09), and dash-dotted (EPS08) lines correspond to 7Pb 
contributions with no shadowing, weak shadowing, moderate 
shadowing, and strong shadowing respectively. Right panel: 
total rapidity distributions (sum of 7p and 7Pb contribu- 
tions) . 



FIG. 10: (Color online) Total rapidity distributions of ex- 
clusive photoproduction of T(ls) in PbPb collisions at the 
LHC in the modified hard sphere density distribution ap- 
proximation. Dashed line depicts result using the MSTW08 
gluon distribution (no nuclear modifications). Solid, dash- 
dotted, and dash-double-dotted lines are results from nuclear- 
modified gluon distributions from EPS09, EPS08, and HKN07 
parton distributions respectively. 



C. Theoretical errors and uncertainties in parton 
distributions 



TABLE VIII: Total cross sections (in fih) for elastic photo- 
production of T(ls) in ultraperipheral PbPb collisions at the 
LHC. 
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well as an efficient discriminator of the different gluon 
shadowing scenarios considered in the current study. 

In order to further exhibit the constraining potentials 
of T photoproduction we display the rapidity distribu- 
tions in Fig. [TU] for the four different gluon shadowing 
scenarios. 

Shadowing remains the relevant nuclear modification 
for practically the entire rapidity range shown in the fig- 
ure, and is markedly manifested in the interval — 2 < y < 
2. Thus rapidity distribution in this interval should be a 
good discriminator of gluon shadowing strength. 



From a consideration of the general structure of the 
quantities (cross sections and rapidity distributions) con- 
sidered in this study, three major sources of theoretical 
errors can be readily identified: 

• accuracy of the relevant expressions for the photon 
flux, 

• higher-order corrections (neglect of terms greater 
than leading order in the case of photoproduction of 
heavy quarks, the accuracy of mopping up higher- 
order corrections with data from HERA in the case 
of production of J /iff and T), 

• uncertainties in the parton distributions. 

Here we give a brief qualitative comment on the third, i.e. 
the effects of uncertainties in parton distributions on our 
calculations. Let us first consider uncertainties in nuclear 
parton distributions (nPDs). This relevant for the dom- 
inant direct contribution in photoproduction of heavy 
quarks and elastic production of vector mesons. Both 
HKN07 HI and EPS09 [H| have facilities for estimating 
error in the determination of a nPDs-dependent quan- 
tity like the cross section due to the uncertainties in the 
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nuclear parton distributions. The HKN07 set achieves 
this through the availability of a set of gradient func- 
tions while EPS09 has thirty error sets to be used in 
conjunction with free proton parton distributions error 
sets. 

For technical reasons we have not computed errors on 
cross sections using either of these two sets. Since this is 
a continuation of our earlier study, we have maintained 
the usage of the MSTW08 set for the free proton parton 
distributions while the EPS09 error evaluation is better 
handled with the CTEQ sets [l^] of nucleon parton distri- 
butions. In the case of HKN07 we have used the nuclear 
modifications in our calculations whereas the gradient 
terms are for the parton distributions proper. We thus 
advocate that the spread in the results involving the nu- 
clear parton distributions be taken as loosely indicative 
of the effect of the uncertainties in the nuclear PDs. 

We now consider the resolved component of photopro- 
duction of heavy quarks where photon parton distribu- 
tions are required. Among the available sets of photon 
parton distributions, only the CJK set has the facility 
to compute errors from uncertainties in photon PDs in 
terms of a tolerance parameter T [6!l[73j], with the CJK2 
uncertainty bands generated with T = 5. Thus, in princi- 
ple, one needs only the CJK2 best fit and the associated 
error sets in order to compute the errors in the resolved 
contribution due to uncertainties in 7PDS. The problem 
here is that the SaSld distributions seems to lie consis- 
tently outside the CJK2 parton distributions error bands 
(further details can be found in [73J). A reflection of this 
trend can be seen in our results for the resolved com- 
ponent where the CJK2 values are consistently about 
a factor of two greater than the SaSld values. As a 
test we have used the CJK2 error sets in conjunction 
with EPS09 to evaluate the uncertainty on the cc re- 
solved cross section in PbPb collisions, with the result 
a (jjK2 = 212.6 ± 10. 7T mb. The corresponding cross 
sections using SaSld and GRV are cr^ sid = 101.1 mb 
and <J t giiv = 169.3 mb respectively. Thus the GRV result 
is already within the error limits of the CJK for T ~ 4 
whereas one needs T ~ 10.4 in order to accommodate the 
SaSld result. In view of this situation and since we have 
no cogent reason to exclude the SaSld distributions, we 
have presented results from all three sets of 7PDS and 
adopt the viewpoint espoused above that the spread in 
results should be taken rather loosely as an indicator of 
the effects of the uncertainties in the 7PDS. 



duction of vector mesons (J/^f and T(ls)) in ultrape- 
ripheral proton-lead (pPb) and lead-lead (PbPb) colli- 
sions at LHC. Both the dominant direct component in 
photoproduction of heavy quarks and elastic production 
of heavy mesons are dependent on nuclear gluon distri- 
butions, and could therefore be potentially useful in con- 
straining modifications such as shadowing and antishad- 
owing in nuclear gluon distributions. The resolved com- 
ponent in photoproduction of heavy quarks depends on 
the distributions of light quarks and gluons in both pho- 
tons and nuclei. Different sets of both nuclear and photon 
parton distributions with different attributes have been 
utilized. 

In photoproduction of heavy quarks the parton dis- 
tribution dependence is linear and different modifica- 
tions are superimposed due to the integration over the 
momentum fraction x. Despite these limitations, both 
cross sections and rapidity distributions for cc in PbPb 
collisions manifest appreciable sensitivity to shadowing 
around midrapidity and a slight sensitivity to antishad- 
owing at more forward and backward rapidities. Thus cc 
photoproduction offers good constraining potential for 
shadowing, and a somewhat less potential for antishad- 
owing. Although photoproduction of bb is less sensitive 
to modifications than cc, the influence of shadowing is 
evident around midrapidity, and it thus offers some con- 
straining ability for shadowing. While both cc and bb 
total photoproduction cross sections and rapidity distri- 
butions in pPb collisions show little sensitivity to nuclear 
modifications, the rapidity asymmetry ratios in some se- 
lect intervals do exhibit significant sensitivity. The re- 
solved components are appreciable, especially for bb and 
are heavily dependent on the choice of 7PDS. Thus it 
seems feasible that they could be of some use in con- 
straining photon parton distributions. 

The outlook for constraining gluon shadowing is even 
better in the case of vector meson production. Here the 
quadratic dependence on gluon modifications makes elas- 
tic photoproduction of vector mesons particularly attrac- 
tive for constraining purposes. The cross sections and 
rapidity distributions for both J and T(ls) photopro- 
duction in PbPb collisions exhibit very good sensitivity 
to gluon shadowing. Thus both offer remarkable poten- 
tial in constraining the shadowing component of nuclear 
gluon distributions. This is also true for T production 
and to a lesser extent J^J production in pPb collisions. 



IV. CONCLUSIONS 

In the present study we have considered photoproduc- 
tion of heavy quarks (cc and bb) and elastic photopro- 
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